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ABSTRACT 

PlWI-interacting RNAs (piRNAs) defend the genome against transposon activity in animal gonads. The Hsp90 chaperone 
machinery has been implicated in the piRNA pathway, but its exact role remains obscure. Here, we examined the effect of 17- 
N-allylamino-17-demethoxygeldanamycin (17-AAG), an Hsp90-specific inhibitor, on the piRNA pathway. In the silkworm 
ovary-derived BmN4 cells, 17-AAG treatment reduced the level of piRNAs and PIWI proteins. In vitro, the 5'-nucleotide 
preference upon precursor piRNA loading was compromised by 1 7-AAG, whereas 3'-end trimming and 2'-0-methylation were 
unaffected. Our data highlight a role of Hsp90 in accurate loading of precursor piRNAs into PIWI proteins. 

Keywords: chaperone; Hsp90; PIWI; PlWI-interacting RNAs; piRNAs 



INTRODUCTION 

PIWI proteins bind to 23-30 nucleotides (nt)-long small 
RNAs called PlWI-interacting RNAs (piRNAs) to form the 
RNA-induced silencing complex (RISC) (Klattenhoff and 
Theurkauf 2008; Ghildiyal and Zamore 2009; Malone and 
Hannon 2009). piRNA- RISC plays an important role in si- 
lencing transposable elements in animal gonads. Indeed, in 
many organisms, mutations in PIWI genes result in de-silenc- 
ing of transposons, causing defects in germline development. 

Unlike small interfering RNAs (siRNAs) and microRNAs 
(miRNAs), piRNAs are born single stranded in a Dicer-inde- 
pendent manner (Vagin et al. 2006; Houwing et al. 2007). In 
the current model for piRNA- RISC formation, single- strand- 
ed primary transcripts are first produced from piRNA-gener- 
ating loci. Then, these transcripts are likely fragmented into 
precursor piRNAs (pre-piRNAs) that are longer than mature 
piRNAs, through multiple steps including endonucleolytic 
cleavage by Zucchini (Brennecke et al. 2007; Klattenhoff 
et al. 2009; Ipsaro et al. 2012; Kawaoka et al. 2012; Muerdter 
et al. 2012; Nishimasu et al. 2012). A subset of PIWI proteins 
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(e.g., Siwi in the silkworm) preferentially binds to pre-piRNAs 
with 5' U (1U) (Brennecke et al. 2007; Gunawardane et al. 
2007; Kawaoka et al. 2009). After PIWI loading, 3' ends of 
pre-piRNAs are trimmed by an enzyme named Trimmer to 
the mature piRNA length (Kawaoka et al. 2011). Trimmer is 
biochemically characterized as a Mg 2+ -dependent 3'-to-5' 
exonuclease, but its identity remains unknown. Finally, cou- 
pled with trimming, 3' ends are 2'-0-methylated by Henl 
(Horwich et al. 2007; Kirino and Mourelatos 2007; Ohara 
et al. 2007; Saito et al. 2007; Kawaoka et al. 201 1), generating 
mature RISC containing primary 1U piRNAs. This is called 
the primary processing pathway. Via RISC, primary piRNAs 
cleave their complementary target RNAs across from posi- 
tions 10 and 11. The 3' fragments of the cleavage products 
are then incorporated into another subset of PIWI proteins 
(e.g., BmAgo3 in the silkworm), and their 3' ends are like- 
ly trimmed to generate secondary piRNAs, which often have 
adenine at position 10 (10A) (Brennecke et al. 2007; 
Gunawardane et al. 2007; Kawaoka et al. 2009, 2011). Such a 
cleavage-dependent piRNA biogenesis is called the ping- 
pong amplification loop or secondary processing pathway 
(Brennecke et al. 2007; Gunawardane et al. 2007). In insects, 
1U primary piRNAs are often antisense to transposons, and 
10A secondary piRNAs are sense to transposons, but some 
transposon-derived piRNAs show opposite strand polarity. 
In theory, 10A secondary piRNAs can further direct produc- 
tion of 1 U secondary piRNAs, although they are indistinguish- 
able from 1U primary piRNAs. 
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Recent reports have implicated Heat shock protein 
90 (Hsp90) and its cochaperones in the piRNA pathway 
(Specchia et al. 2010; Gangaraju et al. 2011; Olivieri et al. 
2012; Preall et al. 2012; Xiol et al. 2012). In flies, disruption 
of Hsp90 causes reduction of a subset of piRNAs and de-si- 
lencing of a number of transposons (Specchia et al. 2010). 
Hsp90 interacts with Piwi, one of the three fly PIWI proteins, 
and regulates the phosphorylation level of Piwi without affect- 
ing its protein stability (Gangaraju et al. 2011). Moreover, an 
Hsp90 cochaperone, Shutdown, plays a critical role in both 
primary and secondary piRNA biogenesis in flies (Olivieri 
et al. 2012; Preall et al. 2012), while Fkbp6, a mouse ortholog 
of Shutdown, is exclusively required for secondary piRNA 
biogenesis (Xiol et al. 2012). In the silkworm ovary-derived 
cell line BmN4, inhibition of Hsp90 causes accumulation of 
16-nt RNAs in BmAgo3 (but not in Siwi), which presumably 
correspond to the 5' fragments of Siwi-directed cleavage prod- 
ucts (Xiol et al. 2012). Given these diverse 
observations, the exact role of Hsp90 in 
the piRNA pathway still remains obscure. 

In this study, we investigated the role Amount loaded: 

of Hsp90 in the piRNA pathway using siwi 
BmN4 cells and a BmN4- cell- derived in 
vitro system that recapitulates a part of 
the primary piRNA biogenesis (Kawaoka 
et al. 2011). Our data highlight a role 
of Hsp90 in accurate loading of pre- 
piRNAs into PIWI proteins. 



proteins was recapitulated by using cells stably expressing 
Flag-Siwi or Flag-BmAgo3 (Fig. 1A). In contrast, quantitative 
PCR (qPCR) analysis revealed that 17-AAG treatment did not 
alter the steady-state mRNA levels of these two genes (Fig. 
IB). Thus, the reduction of Siwi and BmAgo3 occurred at 
the protein level. 

Inhibition of Hsp90 reduces piRNAs in vivo 

To determine if 17-AAG treatment affects piRNA population, 
we analyzed RNAs from BmN4 cells treated with 17-AAG or 
DMSO. The overall piRNA level, simply visualized by SYBR 
Gold staining, was reduced by 17-AAG (Supplemental Fig. 
S3). We then immunoprecipitated Siwi and BmAgo3 and 
visualized the bound RNAs by radiolabeling (Fig. 1C,D). As 
reported previously (Xiol et al. 2012), 16-nt RNA popula- 
tion was detected in the BmAgo3-immunoprecipitates upon 
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RESULTS AND DISCUSSION 

Inhibition of Hsp90 reduces the two 
PIWI proteins in vivo 

To investigate the role of Hsp90 in piRNA 
biogenesis, we first sought to deplete 
Hsp90 in BmN4 cells by RNAi. How- 
ever, the knockdown efficiency was ex- 
tremely low (data not shown), and we 
decided to inhibit the Hsp90 activity 
by 1 7-allylamino- 1 7-demethoxygeldana- 
mycin (17-AAG), a specific Hsp90 antag- 
onist (Kamal et al. 2003). We treated 
BmN4 cells with 17-AAG (dissolved in 
DMSO) for 48 h. In this condition, cell vi- 
ability was comparable to that of the cells 
treated by DMSO only (Supplemental 
Fig. SI). As previously reported (Xiol 
et al. 2012), 17-AAG treatment induced 
de- silencing of a number of transposons, 
especially Pao (Supplemental Fig. S2). 
Moreover, protein levels of both silk- 
worm PIWI proteins, Siwi and BmAgo3, 
were severely decreased by 17-AAG treat- 
ment (Fig. 1A). Reduction of the PIWI 
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FIGURE 1 . Inhibition of Hsp90 reduces Siwi and BmAgo3 proteins in vivo. (A) Wild-type (WT) 
BmN4 cells, Flag-Siwi-expressing cells, and Flag-BmAgo3 expressing cells were treated with DMSO 
or 10 uM 17-AAG for 48 h. Total cell lysates were subjected to Western blot analysis with anti- Siwi, 
anti-BmAgo3, anti-Flag, and anti-Actin antibodies. Actin served as a loading control. (B) The levels 
of Sim and BmAgo3 mRNAs in the BmN4 cells treated with DMSO or 1 0 uM 1 7-AAG for 48 h were 
measured by qPCR. rp49 served as an internal control. The average ± standard deviations from 
three independent trials are shown. (C,D) Wild-type (WT) BmN4 cells, Flag-Siwi, or Flag- 
BmAgo3 stably expressing cells were treated with DMSO or 1 0 uM 1 7-AAG for 48 h. The cell lysates 
were subjected to immunoprecipitation with anti-Flag antibodies and analyzed by Western blot 
(C). RNAs bound to Flag-Siwi or Flag-BmAgo3 were purified from the immunoprecipitates 
and detected by the 3'-end radiolabeling with [a- 32 P]cordycepin 5' -triphosphate (D). Siwi- and 
BmAgo3-bound piRNAs were reduced in proportion to the reduction of the Piwi proteins. (*) 
An ~16-nt short RNA accumulation in BmAgo3 immunoprecipitates in the presence of 17-AAG. 
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17-AAG treatment (Fig. ID). In addition, 
we observed a reduction of Siwi- and 
BmAgo3-bound piRNAs to a degree con- 
sistent with the decrease of these proteins 
by 17-AAG (Fig. 1C,D). To further ana- 
lyze how Hsp90 inhibition affects the 
piRNA abundance and composition, we 
carried out deep sequencing of total small 
RNAs and Siwi- and BmAgo3-bound 
piRNAs with or without 17-AAG treat- 
ment (Supplemental Table 1). In the silk- 
worm piRNA pathway, Siwi preferentially 
incorporates 1U piRNAs, while BmAgo3 
does not show a 5' -nucleotide preference 
(Kawaoka et al. 2009, 2011). Instead, 
BmAgo3-bound secondary piRNAs are 
enriched for 10A. For simplicity, our anal- 
yses excluded piRNAs with both 1U and 
10A, thereby focusing on "1U (but not 
10A)" and "10A (but not 1U)" piRNAs 
that matched the silkworm genome. 
When normalized to 5S ribosomal RNA 
fragments (see Materials and Methods), 
17-AAG treatment decreased 1U piRNAs 
by ~60% and 10A piRNAs by ~40% in 
their abundance in the total small RNA 
library (Fig. 2A), agreeing with the reduc- 
tion of Siwi- and BmAgo3-bound piRNAs 
on the gel ( Fig. 1 D ) . On the other hand, se- 
quencing analysis of Siwi- and BmAgo3- 
bound piRNAs showed no substantial 
change in the overall 1U/10A nucleo- 
tide composition by 17-AAG treatment 
(Fig. 2B). Classification of Siwi- and 
BmAgo3-bound piRNAs by individual 
transposons from which they derived 
also exhibited little overall change in their 
ratio between 1U and 10A piRNAs (Fig. 
2C), although ~35% of Siwi-bound trans- 
poson piRNAs and ~48% of Ago3-bound 
transposon piRNAs pointed to a slight but 
significant (P < 0.01 by Fisher's exact test) 
loosening of their specific nucleotide 
preferences (Fig. 2C). The length distribu- 
tion of Siwi-bound piRNAs (27-28 nt) 
and BmAgo3-bound piRNAs (26-27 nt) 
was also unaffected by 17-AAG (Fig. 
2D). However, in the total piRNA pool, the peak of the 
length distribution showed a shift from 28 nt to 27 nt by the 
17-AAG treatment (Fig. 2E). These results together suggest 
that both Siwi- and BmAgo3 -bound piRNAs were reduced 
by Hsp90 inhibition with little change in their compositions, 
and the reduction was more pronounced for Siwi-bound 
piRNAs. 
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FIGURE 2. Inhibition of Hsp90 reduces Siwi- and BmAgo3-bound piRNAs in vivo. (A) 
Abundance of 1U (but not 10A) and 10A (but not 1U) piRNAs in the total small RNA library 
with DMSO or 17-AAG treatment. Reads were normalized to 5S ribosomal RNA fragments. 1U 
piRNAs and 10A piRNAs were reduced by ~60% and ~40%, respectively, by 17-AAG. (5) The 
1U/10A composition of Siwi- and BmAgo3-bound piRNAs was essentially unaffected by 17-AAG 
treatment. (C) Comparison of the 1U/10A ratio of piRNAs for each transposon with DMSO (x- 
axis) or 17-AAG treatment (y-axis). (Red circles) Siwi- and (blue circles) BmAgo3-bound 
piRNAs; the color density reflects the P-value of a Fisher's exact test for 2 x 2 contingency tables 
consisting of the numbers of 1U and 10A piRNA reads under DMSO and 17-AAG treatment con- 
ditions. The change was small overall, but statistically significant changes of individual transposon- 
derived piRNAs pointed to loosening of the specific 1U and 10A biases of Siwi- and BmAgo3-bound 
piRNAs, respectively. (D) The length distribution of Siwi-bound piRNAs (27-28 nt) and BmAgo3- 
bound piRNAs (26-27 nt) was essentially unaffected by 17-AAG. (E) Length distribution and the 
abundance of piRNAs in the total small RNA library with DMSO or 1 7-AAG treatment. Reads 
were normalized to 5S ribosomal RNA fragments. The shift in the length distribution indicates 
that the reduction of Siwi-bound piRNAs was more pronounced than BmAgo3-bound piRNAs. 



Hsp90 facilitates accurate pre-piRNA loading in vitro 

A previous study demonstrated that unloaded, RNA-free 
Argonaute2 (Ago2) is destabilized in the absence of Hsp90 ac- 
tivity in mammalian cells, whereas Ago2-RISCs "pre-loaded" 
with small RNAs remain highly stable even after Hsp90 inhi- 
bition (Johnston et al. 2010). The overall reduction of Siwi- 
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and BmAgo3-bound piRNAs as well as their steady-state pro- 
tein levels (Figs. 1, 2) is consistent with the view that, upon 17- 
AAG treatment, de novo biogenesis of piRNAs was impaired 
and the resultant RNA-free PIWI proteins were largely de- 
graded, while a preexisting pool of piRNA-RISCs remained 
stable. A part of the piRNA biogenesis pathway — loading of 
pre-piRNAs into Siwi, and trimming and 2'-0-methylation 
of their 3' ends — can be recapitulated in vitro (Kawaoka 
et al. 201 1). Thus, we examined the effect of Hsp90 inhibition 
at each of these steps using this system. Since Siwi-bound 
piRNAs predominantly begin with U in vivo, the 1U bias is 
an important indicator of canonical Siwi-RISC assembly in 
vitro. We performed in vitro pre-piRNA loading assay using 
lysate from BmN4 cells expressing Flag- Siwi and synthetic 
50-nt single-stranded RNAs beginning with U, A, G, or C 
(1U/1A/1G/1C-50 RNAs) (Fig. 3A). We incubated the 
lysate with 5'- 32 P-radiolabeled 1U/1A/1G/1C-50 RNAs in 
the presence or absence of 17-AAG and then immunoprecip- 
itated Flag-Siwi by anti-Flag antibody. The efficiency of Flag- 
Siwi immunoprecipitation was unaffected by 17-AAG (data 
not shown). Finally, we extracted RNAs and measured the 
amount of Siwi-bound RNAs. In the absence of 17-AAG, 



Siwi preferentially incorporated 1U-50 RNAs, recapitulat- 
ing canonical pre-piRNA loading (Fig. 3B,C). In the presence 
of 17-AAG, however, Siwi largely lost the 1U preference 
(Fig. 3B,C). Geldanamycin, another Hsp90 inhibitor, showed 
a similar effect (Supplemental Fig. S4). Thus, Hsp90 is re- 
quired for Siwi to load pre-piRNAs properly. In contrast, 
17-AAG did not affect the 3' -end trimming of Siwi-loaded 
pre-piRNAs or 2'-0-methylation (Fig. 3D,E). Therefore, 
Hsp90 is specifically required for accurate pre-piRNA loading 
into Siwi but not for the subsequent maturation steps (Fig. 
3F). This is consistent with our in vivo observation that the 
abundance of Siwi-bound piRNAs was reduced without af- 
fecting their length (Figs. ID, 2A,D). 

Accumulation of 16-nt RNAs in Ago3 upon 17-AAG treat- 
ment (Fig. ID; Xiol et al. 2012) indicates that Hsp90 functions 
not only in primary loading of pre-piRNAs into Siwi but 
also in the ping-pong amplification of secondary piRNAs. 
In the current model of the piRNA biogenesis, BmAgo3 should 
not autonomously incorporate single-stranded RNAs, but 
should rather exclusively receive the cleavage products by 
Siwi. Unfortunately, there has been no in vitro system that re- 
capitulates such an RNA hand-over mechanism. However, 
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FIGURE 3. Hsp90 facilitates accurate pre-piRNA loading in vitro. (A) Experimental procedure for B. (B) Fifty-nucleotide 5'-radiolabeled RNAs with 
5' U, A, G, or C were incubated in lysates prepared from wild-type BmN4 cells (naive) or Flag-Siwi- or Flag-BmAgo3-expressing cells in the presence 
or absence of 1 mM 17-AAG. After immunoprecipitation with anti-Flag antibody, Siwi- or BmAgo3-bound RNAs were detected. (C) Quantification 
for B. The level of Siwi- or BmAgo3-bound 1U-50 RNA in DMSO control was normalized to 1. The average ± standard deviations from three inde- 
pendent trials are shown. 17-AAG abolished the 1U preference of Siwi, while promoting promiscuous incorporation into BmAgo3. (D) Experimental 
procedure for E. (E) 1U-50 RNA loaded into Flag-Siwi was subjected to the trimming assay in the presence or absence of 1 mM 17-AAG. The 27- to 
28-nt trimming products were generated even in the presence of 17-AAG and were resistant to NaI0 4 -mediated oxidization followed by (3-elimination, 
which is indicative of 2'-0-methylation. As a control, addition of 1 mM S-adenosyl homocysteine (SAH) inhibited 2'-0-methylation of the trimming 
product, which resulted in its faster gel migration after NaI0 4 /(3-elimination. (F) A function of Hsp90 in the primary piRNA pathway. Hsp90 is re- 
quired for accurate loading of pre-piRNAs into Siwi, but not for subsequent 3'-end trimming and 2'-0-methylation. 
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in vitro, BmAgo3 can incorporate 50-nt RNAs to some ex- 
tent, without showing significant 5' -nucleotide preference 
(Kawaoka et al. 201 1), which likely represents a marginal ac- 
tivity of BmAgo3 to load primary piRNAs. Intriguingly, 17- 
AAG enhanced the overall incorporation of 50-nt RNAs 
into BmAgo3 (Fig. 3B,C). The loss of the 1U bias of Siwi 
and the enhanced overall RNA incorporation of BmAgo3 
by 17-AAG treatment in vitro may suggest that Hsp90 func- 
tions to ensure the distinct RNA-loading mechanisms of Siwi 
and BmAgo3. Supporting this idea, both the 1U bias of Siwi- 
bound piRNAs and the 10A bias of BmAgo3-bound piRNAs 
were slightly but significantly loosened when BmN4 cells 
were treated with 17-AAG (Fig. 2C). Notably, in fly ovaries, 
depletion of Aubergine, of which Siwi is the silkworm coun- 
terpart, caused accumulation of cytoplasmic foci containing 
unloaded Ago3, Krimper (a Tudor-domain protein), 
Hsp90, and Shutdown (Malone et al. 2009; Olivieri et al. 
2012). Similarly, in the fly ovarian somatic cells, where only 
the primary processing pathway operates, ectopically ex- 
pressed fly Ago3 failed to incorporate piRNAs, and such un- 
loaded Ago3 and Shutdown accumulated in Krimper- 
enriched cytoplasmic foci called Krimp bodies (Olivieri 
et al. 2012). Thus, the Hsp90 chaperone machinery together 
with Krimper may sequester unloaded Ago3 protein from 
promiscuously incorporating primary piRNAs by maintain- 
ing its proper structure, until Ago3 is docked on the ping- 
pong platform. Further investigation is warranted to address 
these hypotheses. 

MATERIALS AND METHODS 
Cell lines 

BmN4 cells were cultured at 27°C in IPL-41 medium (Applichem) 
supplemented with 10% fetal bovine serum. 

General methods for in vitro experiments 

Preparation of 40x reaction mix (containing ATP, ATP regeneration 
system, and RNase inhibitor), lysis buffer (30 mM HEPES-KOH at 
pH 7.4, 100 mM KOAc, 2 mM Mg(OAc) 2 ), and 2x PK buffer has 
been previously described in detail (Haley et al. 2003). 

In vitro assay for piRNA biogenesis 

In vitro pre-piRNA loading, trimming, NaI0 4 -mediated oxidation, 
and P-elimination were performed as described previously 
(Kawaoka et al. 2011). 17-AAG (Sigma-Aldrich) was used at a final 
concentration of 1 mM. 

Western blotting 

Western blotting was performed as described previously with anti- 
Flag (1:1000; Sigma-Aldrich), Actin (1:3000; Santa Cruz), Siwi 
(1:3000), and BmAgo3 (1:3000) antibodies (Kawaoka et al. 2009). 



RNA extraction, reverse transcription, quantitative 
PCR, and 3 -end labeling 

Total RNAs were extracted by using TRIzol reagent (Invitrogen) or 
the miRVana miRNA isolation kit (Ambion) according to the man- 
ufacturers' instructions. Reverse transcription for mRNAs and the 
following quantitative PCR analyses were performed essentially as 
described previously (Kawaoka et al. 2008). The primer sequences 
for quantitative PCR analyses are listed in Supplemental Table 
2. RNA 3'-end labeling was performed with 2 uM [a- 32 P]cordycepin 
5'-triphosphate (Perkinelmer) and yeast poly(A) polymerase (USB) 
for 30 min at 37°C. 

piRNA library construction and bioinformatic analyses 

piRNA libraries were constructed as described previously (Kawaoka 
et al. 2009). To analyze the change of piRNA abundance by 17-AAG, 
we used the number of reads matching 5S ribosomal RNA fragments 
for normalization because they were relatively unsusceptible to 17- 
AAG treatment; miRNAs were unsuitable for normalization pur- 
pose, since Hsp90 plays a critical role in their biogenesis (Iwasaki 
et al. 2010; Johnston et al. 2010). When analyzing the change of 
the nucleotide composition of Siwi- and BmAgo3-bound piRNAs, 
the total genome-mapping reads were used for normalization. 
Bioinformatic analyses were performed with the aid of in-house 
UNIX and R programs, which will be provided upon request. 

DATA DEPOSITION 

The piRNA deep sequencing data used in this study are available 
under the accession number of DRA000943 (DDBJ). 

SUPPLEMENTAL MATERIAL 

Supplemental material is available for this article. 
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